Response of a particle in a one-dimensional lattice to an applied force: Dynamics of 

the effective mass 
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We study the behaviour of the expectation value of the acceleration of a particle in a one- 
dimensional periodic potential when an external homogeneous force is suddenly applied. The theory 
is formulated in terms of modified Bloch states that include the interband mixing induced by the 
force. This approach allows us to understand the behaviour of the wavepacket, which responds with 
a mass that is initially the bare mass, and subsequently oscillates around the value predicted by 
the effective mass. If Zener tunneling can be neglected, the expression obtained for the acceleration 
of the particle is valid over timescales of the order of a Bloch oscillation, which are of interest for 
experiments with cold atoms in optical lattices. We discuss how these oscillations can be tuned in 
an optical lattice for experimental detection. 

PACS numbers: 37.10.Jk, 37.10.Vz, 72.20.-i 



I. INTRODUCTION 

In the absence of scattering due to impurities and 
phonons, the wavepacket associated with a crystal elec- 
tron accelerates in response to an external force F (ho- 
mogeneous in space and constant in time) as a particle 
with an effective mass. This observation is justified by 
the effective mass theorem [1], which in the simple case 
of a one-dimensional lattice is 



F 



(1) 



where (a) is the expectation value of the acceleration of 
the wavepacket and m*(fc) is inversely proportional to 
the curvature of the band energy, 



h 2 dk 2 



(2) 



However, Pfirsch and Spenke [2] argued that (a) does 
not satisfy (1) at all times when the force is suddenly 
applied; instead, the initial response of the expectation 
value of the acceleration is characterized by the bare 
mass of the electron. There should be subsequent oscilla- 
tions around the value calculated from the usual effective 
mass (2), which have a period inversely proportional to 
the energy gap. They die off after a characteristic time 
that roughly decreases with decreasing bare mass and 
lattice constant [2]. 

In typical solid-state systems, the femtosecond scale 
of this characteristic time and the scattering due to 
impurities and phonons make it difficult to observe 
these oscillations. Atoms in optical lattices constitute a 
much simpler system, where the timescale is expected 
to be longer and decoherence can be minimized. In 



these systems a constant force can be introduced by 
accelerating the lattice uniformly so that the atoms 
experience an inertial force in the lattice frame. Many 
experiments with ultracold atoms have been carried out 
to examine interesting phenomena caused by the band 
structure in one-dimensional optical lattices [3-8]. Here 
we discuss the possibility of detecting the oscillatory 
behaviour of the expectation value of the acceleration 
in ultracold atoms in a one-dimensional optical lattice 
when an inertial force is suddenly applied. 

Assuming that Zener tunneling is not significant, 
we write an approximate semianalytical expression 
for the expectation value of the acceleration valid for 
times as long as one Bloch period (see expression (48)). 
Our approach is based on constructing a wavepacket 
from modifed Bloch states that take into account the 
first order interband mixing due to the external force. 
These modified Bloch states are constructed using a 
diagonalization scheme to decouple the bands to some 
given order in the external force [9], following ideas first 
introduced by Adams [10-12], Kane [13] and Wannier 
[14]. To first order, we write an explicit expression 
for the desired expectation value of the wavepacket 
acceleration. We find that at early times it reduces 
to the result found by Iafrate and Krieger in their 
study of the motion of crystal electrons shortly after a 
dc field is applied [15-18]. We discuss the validity of 
our perturbation approach, comparing with the results 
obtained from a full numerical calculation, and confirm 
that our approach is still valid for longer times. We 
identify the effect of changing the different physical 
parameters of the optical lattice, such as the bare 
mass of the atoms and the lattice constant, and find 
some sets of parameters that we believe would lead to 
experimentally measurable effects. 
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To the best of our knowledge, the type of oscillations 
described here have not been measured experimentally. 
Although the observation of this phenomenon in optical 
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lattices is an interesting problem on its own, we believe 
it can also help to understand how such oscillations 
could be detected in solid-state systems where the 
femtosecond timescales are becoming more accessible 
through ultrafast pulses [19]. 

In Section II we summarize the strategy to calculate 
the expectation value of the acceleration, in terms of 
modified Bloch states, and illustrate its behaviour for a 
one-dimensional Mathieu potential with parameters ad- 
justed to resemble those of the band structure of a semi- 
conductor. The method used for the full numerical calcu- 
tion is sketched in Section III. In Section IV, we show an 
application of the formalism to cold atoms in an optical 
lattice with several examples, illustrating the relevant pa- 
rameters that control the oscillations. We present some 
conclusions in Section V. 



II. THEORETICAL FRAMEWORK 

We consider a particle in an infinite one-dimensional 
lattice described by a wavepacket with spread in quasi- 
momentum smaller than the extension of the Brillouin 
zone. In this section we describe the formalism employed 
to include the effect of an external homogeneous force 
acting on the particle. 



A. Hamiltonian and modified Bloch states 

The Hamiltonian that describes the system is 

H = Ho- F(t)x, (3) 

where H corresponds to the unperturbed part of the 
Hamiltonian, 



(4) 



The potential V(x) has period b. We will treat —F{t)x 
as a perturbation, but in a particular sense that will be 
described later. 

If we wanted to solve this problem in the crystal mo- 
mentum representation we would expand the wavefunc- 
tion of the system ^(x, t) as a wavepacket of Bloch func- 
tions, -0„(fc,x), 



= X] dk c n (k,t)i[> n (k,x), 



(5) 



where the integral is over the Brillouin zone. The Bloch 
functions can be written as 



ip n (k,x) = -=u n {k,x)e ikx , 

V Z7T 



(6) 



where u n (k, x) has the periodicity of the one-dimensional 
lattice. The Bloch functions are eigenstates of the un- 
perturbed Hamiltonian W with energies S n (k). They do 
not diagonalize the full Hamiltonian when we include the 
force term, which destroys the periodicity of the system 
and is divergent for \x\ — >■ oo. Therefore, some care is 
required to express the position operator in the crystal 
momentum representation [20] 

/ V^O, x)xip n > (k', x)dx = S nn > I -i—6(k ~ k ') 

+ 6(k-k% ln ,(k), (7) 
where t; nn '(k) is the Lax connection [21] 



f b d 

Un'(k) = J u* n (k,x)i—u n '(k, 



. dx 
x) T . 



(8) 



In the crystal momentum representation the Hamiltonian 
(3) becomes 

/oo 
iP* n (k,x)Hip n ,(k',x)dx = H nn ,{k;t)S{k-k')- 
-OO 

Snn'F(t)-(-i-^S(k-k')j, (9) 

where we have introduced the matrix elements 

H nn ,{k-t) = £*(k;t)5 nn , - F(t) ■ Can-(k). (10) 

Here £^(k;t) denotes the band energy renormalized by 
the diagonal part of the Lax connection, 



£*(k;t) = E n {k)-F(t)-S nn (k). 



(11) 



The time evolution of the amplitudes c n (k, t) is given by 
the time-dependent Schrodinger equation 

d 

ih—c n (k,t) = y~] H nn > (k; t)c n > (k, t)- 



iF(t)- —c n (k,t). (12) 



Now we look at the case where the force F(t) is con- 
stant for all times, that is 



F(t) = F, for t E (-00,00). 



(13) 



Consider a general unitary transformation U n > n (k) of the 
Bloch states, tjj n (k,x), 



i(k,x) = ^2ip n ,(k,x)U n , n (k). 



(14) 



In terms of these modified Bloch states, the wavepacket 
(5) becomes 

*(M) = J2 I dkb n (k,t)<j> n (k,x), (15) 
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where 



where 



b n {k,t) = Y,K-n{k)c nl {k,t). (16) 

n' 

From (12) the evolution of b n (k,t) is given by 
ih^b n (k,t) = E ( E U* mn {k)H mm ,{k)U m , n ,{k)- 

n' mm' 
m ' 

(17) 

In order to decouple the amplitude for the band n from 
the rest of the bands, it would be desirable to find a trans- 
formation U n ' n (k) that diagonalizcs the term in paren- 
theses in (17), 

^ ' U mn {k)H mTn i(k)U m ' n >(k)— 

mm' 

E U *mn{k)iF ■ ^U mn ,{k) = S nn ,W n (k). (18) 



If this equation could be solved, the solution to (17) 
would be 

b n (k, t) = b n (k - ±Ft, 0)e-* /o w ^ k ~ F (t-t-)/h)dt>_ (ig) 



Wannier proved that it is possible to find a transfor- 
mation that satisfies (18) using an expansion of (14) in 
powers of F, and described a recurrence procedure to 
construct it [14]. Such a power series expansion is only 
valid if Zener tunneling between bands is not significant 
[9, 14]. This is appropriate for the physical situations we 
will discuss here, where the force is small enough so that 
the wavepacket remains essentially in one band. In the 
Appendix we summarize Wannier's procedure. To first 
order in F, the unitary transformation U n i n (k) is 



ln {n,t) = l / £*(K+\Ft')dt' (24) 



U n 'n{k) ~ Sn'n "f~ A n ' n (Jt), 



(20) 



where A n > n (k) is a fc-dependent dimensionless parame- 
ter that compares the interband position matrix element 
times the force with the energy difference between bands, 

£-n'n{k) = Sn'(k) — £ n (k), 

A n , n (k)= F ^'^\ 1 - Sn , n ). (21) 

The first order approximation for W n (k) is 

W n (k)n£*(k). (22) 
Hence, the solution (19) valid to first order is [14] 

MM) « b n (k - ^Ft,0)e-^ k -i Ft '*\ (23) 



with k = k — Ft/h. 

Adams and Argyres [10, 11] used the modified Bloch 
states (14) with the first order approximation (20), 

<p n (k,x) w tp n (k,x) +'Y]ip nl (k,x)A n i n (k), (25) 

n' 

to construct states where the expectation value of the 
acceleration, 

d 2 r°° 

(a(t)) = —j **(x,t)x*(x,t)dx, (26) 

obeys the effective mass theorem (1). The parameter 
(21), which controls the mixing between the bands, 
is assumed to be small but it is necessary to dress 
the particle in the periodic potential to establish the 
effective mass. 

The expectation value of the acceleration can be writ- 
ten as 

(a(t)) = ^ + J- f° [p,H ] *(x,t)dx, 

m inm J_ oc 

(27) 

where p denotes the momentum operator and the force 
F(t) can have any time dependence. For a wavepacket of 
the form 

*'jv(M) = J dkb N (k,t)cp N (k,x) (28) 

and a constant force (13), the expectation value (27) be- 
comes 

(a(t)) = ^+f dk\b N (k,t)\ 2 ^-F NN (k), (29) 
where we defined 
■F nri (k) — 

Pmm' 

(k)£ mm >{k)U mW (k) (30) 

mm' 

in terms of the momentum matrix elements 

p nn '{k) = ^- J tpl(k,x)^-^ip n '(k,x)dx. (31) 
To first order in A„/„(fc), F nn i{k) is 



■F nn'{k) ~ — ( £nn' {k)pnn' (k)-\- 



y ' Pnm(k)£nm(k) A mn r (fc) 
m 

^ ] A nm {k)p mn i (k)£ mn i (fc) j. (32) 
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For n = n' = N, equation (32) simplifies to 



•Fjvjv(fc) « F ■ 



m* N (k) 



(33) 



as a result of the well-known sum rule [21] 

m = _2_ Pnn'(k)Pn'n{k) , - 

m*(fc) + m 2^ £ nn ,{k) ' 1 J 

Therefore, the expectation value of the acceleration for 
the wavepacket (28) satisfies 



a(t)) = j dk\b N (k,t)\' 



F(t) 
m* N (k)' 



(35) 



with (f>N(k,x) in the approximation (25). This corre- 
sponds to a particle accelerating with the usual effective 
mass (2) at all times [10, 11]. 



B. Expectation value of the acceleration and 
effective mass 



We are interested in the situation where the 

wavepacket is prepared initially in one band only (de- 
noted by N), 

*(s,0) = J f N (k)tp N (k,x)dk, (36) 

and the external (constant) force is suddenly applied at 
t = 0, that is 



F(t) = Q(t)F, 



(37) 



where O(i) is the Heaviside function. The function /jv(fc) 
is assumed to have a spread smaller than the size of the 
Brillouin zone; for instance, in the numerical calculations 
presented in the next sections we will assume a Gaussian 
shape (centered at k — 0) 



f N (k) 



characterized by the width a. 



1 

e 



(38) 



For the initial wavepacket (36) the term involving the 
commutator in (27) vanishes because the momentum ma- 
trix elements in the crystal momentum representation are 
diagonal in k, 



ll>*(k,x)pipn'(k',x)dx — p nn >(k)5(k - fc'). (39) 



Hence, the particle accelerates initially with its bare 
mass, m [2]. Naively, we might think that the wavepacket 
for t > would be of the form 

$ N (x,t) = J dkc N (k,t)ip N (k,x), (40) 



with cjv(fc,0) = /jv(fc) to satisfy the initial condition 
(36), but then the wavepacket (40) would describe a 
particle accelerating at all times with the bare mass. 
This seems to contradict the effective mass theorem (1); 
however, the contradiction is only apparent because the 
Bloch states of one band only are not the appropriate 
states to describe the particle in the presence of a 
homogeneous force. 

The modified Bloch states (25) include the effect 
of the interband mixing due to the external force (to 
first order) and allow us to construct a wavepacket 
that describes a particle accelerating according to the 
effective mass theorem for a constant force as shown 
in (35). This property suggest that such states form 
an appropriate basis even when the force is turned on 
instantaneously at t = 0; thus, instead of a wavepacket 
of the form (28), we use an expansion of the form (15), 
including the amplitudes b n (k,t) for N, and impose 
the initial condition (36), which forces the particle to 
respond initially with its bare mass. For later times 
we use the result (23) to describe the evolution of the 
amplitudes b n (k,t) in time. 

We find that the initial amplitudes b n (k, t — 0), correct 
to first order in A„/„(fc), are 



Mfc,0)« f N (k) 



(41) 



and 



MM) « -f N (k)A nN (k), for n^N. (42) 
For later times we have 

MM) « M/Oe-^'O (43) 

and 

MM) « -/w( K )A„iv(K)e- l7 " (K '* ) , for n + N, (44) 

where k = k — Ft/h. Within this approximation, the 
amplitude for the modified Bloch state N is of zeroth 
order while the amplitudes for the remaining modified 
Bloch states are of first order. The same is true for the 
amplitudes (correct to first order) for the usual Bloch 
states in the expansion (5), 



cjv(M) « /4«)e" 7K(M) 



(45) 



and 



MM) « f N (K)(A nN (k)e~^^~ 

A nN (K)e-^ K A, for n + N. (46) 



These solutions do not include Zener tunneling from the 
initial band N to the neighbouring bands because the 
amplitudes in the latter remain at higher order in the 
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FIG. 1. First four energy bands (n = 0, 1, 2, 3) for a potential 
strength s = 10. (Color online.) 



expansion [9, 14]. 

Compared to (29), there are some additional terms in 
the new expression for the expectation value of the ac- 
celeration (27). For t > 0, the acceleration is 

( (t)> = [ dk\b n (k,t)\ 2 -T nn (k)+ 

7JI L — ' / 777 

n 

V [ dk-<Rz[b* n {k,t)b n ,(k,t)T nn ,(k)], (47) 

n,n 

where 9te[-] denotes the real part. In the first order ap- 
proximation we find that (47) reduces to 



where Vq is a constant that determines the strength of 
the potential, and fc^ = ir/b identifies half of the Bril- 
louin zone. The eigenvalue problem for the unperturbed 
Hamiltonian H can be rewritten as Mathieu's differen- 
tial equation. The two linearly independent solutions 
of Mathieu's problem can be combined appropriately 
to construct Bloch-type solutions (6) analytical and 
periodic in fc over the Brillouin zone [— fci,fci] for each 
band index [22-25]. 

The energy scale can be characterized by the kinetic 
energy associated with the wavevector fc^, 



Er 



h 2 k 2 L 

2m 



(50) 



and we can write Vq — sEr, where s is a dimensionless 
constant. An example of the band structure calculated 
for s = 10 is shown in Figure 1. For future reference, 
the bands will be labeled in order of increasing energy 
starting with n = for the lowest band. 

Because of the inversion symmetry of the potential 
V(x) and the procedure used to find the Bloch functions, 
the diagonal part of the Lax connection appearing in the 
renormalized band energies (11) vanishes [23] and the 
momentum matrix elements (31) are purely real. For 
such a potential, (48) simplifies to 



F 

m 

£nN(k) 



(a(t)) 

-y 



dk\f N (K)f 



m* N (k) 



TPNn(k)pnN(K)cOS~f Nn (K,t)). (51) 



+ 



_2 £nN(k) 



PNn{k)pnN(K)e 



(48) 



where we have used (32), (33), (43) and (44). As be- 
fore we use k = k — Ft/h and additionally we introduce 
1fNn(K,t) = 7jv(tj t) — Jn(K,t). The first term in (48) 
describes the acceleration with the usual effective mass 
while the second term contains the oscillations predicted 
by Pfirsch and Spenke [2] . Using the sum-rule (34) it can 
verified that (a(t = 0)) = F/m as expected [15, 17, 18]. 



C. Example 

Having presented the general formalism for a one- 
dimensional periodic potential, we will now assume a 
Mathieu potential to illustrate the behaviour of the ac- 
celeration (48). Therefore, we set 



For times short enough after the force is applied, the 
variations of the momentum matrix elements and the 
energy differences as k changes with time can be ne- 
glected and we recover the expression found by Iafrate 
and Krieger for an electron in a dc electric field, using 
a vector potential gauge and accelerated Bloch states 
[15, 17, 18]. However, as will be seen in some examples in 
Section IV, the expression (51) gives a good approxima- 
tion to the behaviour of the acceleration for longer times, 
over which the variations of the momentum matrix ele- 
ments and the energy differences cannot be neglected, as 
long as the wavepacket remains mainly in the original 
band. Such long times might not be of practical inter- 
est for typical solid-state systems where it is difficult to 
maintain coherence over times comparable to the Bloch 
period, 



t b = h/bF, 



(52) 



V{x) = V sm 2 {k L x), 



(49) 



which is the time necessary for a wavepacket to return 
to its original position in the Brillouin zone under the 
action of a constant force [1]. Experiments in optical 
lattices, on the other hand, can access the dynamics 
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FIG. 2. Time dependent effective mass (53) for an electron 
wavepacket (a m 0.2fct) initially at the center of the Bril- 
louin zone in the band N = 2 of a one-dimensional poten- 
tial (49) with s — 10 and b — 0.5 nm, calculated from the 
first order approximation (red solid line); the blue dotted line 
corresponds to the usual effective mass associated with the 
wavepacket. For reference, the time is given both in fem- 
toseconds and in units of the Bloch period (here tb ~ 485 fs). 
The first order approximation is essentially indistinguishable 
from the full numerical calculation result (see Section III). 
We show both results for a smaller time interval in the inset, 
where the red solid line corresponds to the first order approxi- 
mation and the green crosses correspond to the full numerical 
result. The units in the inset are the same as in the main 
plot. (Color online.) 



FIG. 3. Expectation value of the acceleration for the same 
electron wavepacket described in Figure 2 but starting in the 
lowest band, N = 0. The first order approximation (51) cor- 
responds to the red solid line and the value predicted by the 
usual effective mass corresponds to the blue dotted line. The 
light blue slash-dotted line corresponds to the envelope func- 
tion in (58), adjusted to give the right initial value of the 
acceleration. As for the example shown in Figure 2, there 
is excellent agreement between the first order approximation 
and the full numerical calculation result (see Section III). We 
show both results for a smaller time interval in the inset, 
where the red solid line correspond to the first order approxi- 
mation and the green crosses correspond to the full numerical 
result. The units in the inset are the same as in the main 
plot. (Color online.) 



on timescales of the order of a Bloch period relatively 
easily [3, 4], motivating the study of the oscillations of 
the effective mass as the atoms perform a full Bloch 
oscillation. This will be discussed in Section IV. 

We now illustrate some of the features contained in the 
result (51) for an electron in a one-dimensional potential 
(49) with parameters adjusted to resemble those of the 
band structure of a semiconductor; specifically, we choose 
s = 10 and b = 0.5 nm. For this example the initial 
wavepacket is centered at k = in the band N = 2, which 
resembles the first conduction band of a semiconductor 
(see Figure 1). We use a force that corresponds to a 
strong dc field of 1.7x 10 7 V/m and the band gap £21 = 
0) is of the order of an electron volt (here En = 1.5 eV). 
Using the result (51) with these parameters, we plot the 
effective mass of the wavepacket, defined as 



m*(t) 



(53) 



in Figure 2; note that this is to be distinguished from 
the effective mass for an energy band (2). In this case 
we expect the first order approximation to be valid close 
to the center of the Brillouin zone. 



with the bare mass and, afterwards, the effective mass 
(53) oscillates around the usual effective mass (2) aver- 
aged over the extension of the wavepacket. Although 
there is a superposition of various oscillations in the sec- 
ond term of (51), 



2F 
m 



2 E 

£nN(k) 



dk\f N ( K )\ 2 x 

2 PNn(k)pnN(n) COS 7iVn («,*)) (54) 



associated with the different neighbouring bands, Fig- 
ure 2 shows a distinctive oscillation with period 



:(*) 



|£jvri(fc)| ' 



(55) 



given by the energy difference between the initial band, 
N = 2, and its next neighbour band, n = 1, in the region 
near the center of the Brillouin zone. For the example in 
Figure 2 the period is r OS c(0) ~ 1.30 fs. The remaining 
oscillating terms have small contributions because of the 
much larger energy difference with respect to the band 
N = 2. 



This example displays the features discussed by Pfirsch 
and Spenke [2] . Notice that the system behaves initially 



For the same system but with the electron starting 
in the lowest band, N = 0. the acceleration is also 
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described very well by the first order approximation 
(51) as shown in Figure 3. Compared to the previous 
example, the effective mass of the initial band TV = in 
this case is larger because this band is flatter than the 
band N = 2 discussed before. Due to the larger effective 
mass, the average value of the acceleration is smaller 
and the expectation value of the acceleration goes to 
zero several times. Therefore, it is more convenient 
to plot the acceleration instead of the effective mass (53). 

The oscillations decay as a consequence of the spread 
of the wavepacket in k [2] . It is possible to derive a simple 
approximate expression for the envelope of the decaying 
oscillations of the acceleration while the wavepacket is 
moving close to the center of the Brillouin zone. If we 
consider only the contribution from the next neighbour 
band, n, and neglect the change of the momentum matrix 
elements and energy differences in equation (54) [2] we 
can write the approximate expression 



time it takes the oscillations in (58) to reduce by half, 



2F( PnN {0)) 2 
m 2 £ fiN (0) 



dk\f N {k)\ 2 cos ( -£Nn{k)t 



(56) 

The energy difference associated with the frequency of 
the oscillations can be estimated approximating the band 
energies close to k = by parabolas, characterized by the 
usual effective mass, 



£jvn(fc) 



h 2 k 2 



1 



i 



tr^(O) m*(0) 
= ^P) +W0) ' 



+ WO) 



(57) 



where we introduced the reduced effective mass m^(fc) 
associated with the original band N and its closest neigh- 
bour n. Using this approximation and the Gaussian 
quasimomentum distribution (38), equation (56) can be 
evaluated analytically, yielding an oscillating term with 
an envelope function that controls the decay. The ap- 
proximation for (54) reduces to 



^osc(^) 



2F(p fiAr (0)) 2 cos 



fiVn(0)t 



+ i tan" 



1 / Ha A t 



m 2 £n N {0) 



"I 1/4 



(58) 

The oscillations are characterized by a frequency 
|£ivn(0)|//i, but there is an additional contribution to 
the phase of the cosine in (58). This contribution, 
however, is small in the time ranges where the approxi- 
mation is appropriate. 

Naturally, the approximation (58) is limited because 
we have not included the contributions from the remain- 
ing bands or the motion of the wavepacket through the 
Brillouin zone, but it describes appropriately the form of 
the decay as shown, for instance, in Figure 3. We use the 



_ ^|n<(0)| 

Tdccay - ^ 



(59) 



as an estimate for the decay time of the oscillations of 
the acceleration around the value predicted by the usual 
effective mass. For the examples shown in Figures 2 and 3 
the estimated decay times (59) are 1.89 fs and 4.78 fs, 
respectively, which correspond to a small fraction of the 
Bloch period. 



III. FULL NUMERICAL CALCULATION 

According to the discussion in Section II A, the usual 
effective mass can be understood as the result of a 
"dressing" process of the wavepacket in one band with 
small amplitudes over the neighbouring bands [10, 11]; 
the initial response of the wavepacket, according to 
the bare mass (instead of the usual effective mass), 
produces the oscillatory behaviour shown in expression 
(51). However, this simple picture is valid only if the 
wavepacket remains mainly in one band. The approx- 
imate solution (46) for the amplitudes associated with 
the neighbouring bands predicts a small second order 
correction for the probability of finding the wavepacket 
in those bands. Hence, situations where the population 
of the neighbouring bands becomes important due to 
Zener tunneling cannot be described within this scheme. 
In order to verify that the effects of the population of 
neighbouring bands beyond our approximate solution 
can be neglected in our examples, we compare the 
results from (51) with a full numerical calculation that 
solves the time-dependent Schrodinger equation for the 
Hamiltonian (3) assuming that the force has the form 
(37) and the initial condition is given by (36). 

For a full numerical calculation we use the split-step 
operator method in its original version, where the evolu- 
tion due to the kinetic energy term of the Hamiltonian 
is done in Fourier space [26, 27]. The full Hamiltonian, 
given by (3) in the position representation, corresponds 
to the operator H = p 2 /2m, + V(x) — Fx, where p and x 
are the momentum and position operators, respectively. 
It is convenient [4] to transform it to a new version 



H'{t) = S\t)HS{t)-ihS\t) 



according to the unitary operator 



dS(t) 
dt ' 



S(t) = e 



_ p i a (t)p/h-if3(t)x/hi~f(t)/h 



(60) 



(61) 



where a(t) = atf /2, j3{t) = ma^t and j(t) = majif /3. 
The constant is set by the force, according to F = 
—ma,L- The new Hamiltonian H'(t) is time dependent 
for t > and has the same periodicity as the unper- 
turbed Hamiltonian. It can be written in the position 
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FIG. 4. Expectation value of the acceleration of a wavepacket 
(a — 0.2A;l) initially at the center of the Brillouin zone in 
the band N = for a rubidium atom in an optical lattice 
with s = 7 and b = 390 nm, calculated from the first or- 
der approximation (red solid line), the full numerical solution 
(green dashed line) and the usual effective mass (blue dot- 
ted line). The acceleration of the lattice is 24.2 m/s 2 , which 
corresponds to F ~ 0.173 (see (66)). Notice how well the ap- 
proximate calculation reproduces the full numerical solution. 
(Color online.) 



FIG. 5. Expectation value of the velocity of rubidium atoms 
in the same situation described in Figure 4, calculated from 
the first order approximation (red solid line), the full numeri- 
cal solution (green crosses) and the usual effective mass (blue 
dotted line). The inset shows the difference between the first 
order approximation and the usual effective mass prediction. 
The units in the inset are the same as in the main plot. (Color 
online.) 



behaviour of the acceleration (51). 



representation as 

h 2 d 2 



H'(t) 



a + sE R sm 2 (k L {x- -a L t 2 )). (62) 



In the case of an electron in a dc electric field, this trans- 
formation is equivalent to a gauge transformation where 
the new Hamiltonian is written in terms of the vec- 
tor potential, as was done by Krieger and Iafrate [16-18]. 



In experiments with cold atoms in optical lattices, a 
constant and homogeneous force can be introduced by 
uniformly accelerating the lattice with a linear increase 
in time of the frequency difference between the two 
interfering laser fields that create the lattice [3]. The 
appropriate Hamiltonian to describe the system in 
the laboratory frame is H lab (t) = H'(t), introduced 
for convenience in Section III, with now being the 
acceleration of the lattice [41 . 



The results presented in Figures 2 and 3, calculated 
from expression (51), are in almost perfect agreement 
with the full numerical calculation, confirming that (51) 
is appropriate for short times after the force has been 
suddenly applied. In the next section we will show that 
(51) is also an useful estimate for longer times of the 
order of a full Bloch oscillation, which are relevant for 
experiments with cold atoms in optical lattices. 



IV. COLD ATOMS IN OPTICAL LATTICES 

As illustrated in Section II C, the oscillations dis- 
played in Figures 2 and 3 decay in a few femtoseconds 
for electron motion in a semiconductor, a timescale that 
is difficult to resolve experimentally. We will see that 
the corresponding timescales in optical lattices are of 
the order of microseconds. Considering the additional 
advantages of tunability and low decoherence that 
experiments in optical lattices offer, we believe that such 
systems are excellent candidates to study the oscillatory 



For optical lattices the physical scales change dras- 
tically compared to the electron case described in Sec- 
tion IIC, but the behaviour of the oscillations in (51) 
is qualitatively the same. In order to compare the two 
cases, it is appropriate to consider the ratio of the initial 
oscillation period (55) to the Bloch period, 



Tosc(O) bF 



ttF 



(63) 



T B \£Nn(0)\ |£ Wfl (0)|' 

and the ratio of the decay time (59) to the Bloch period, 
Tdccay _ >/l5 6F|m?&(0)| _ >/l5 F|m5&(0)| 



TB 



(64) 



In (63) and (64) we included the ratios in terms of the 
scaled wavevector, energy and mass: 

fc = A^(fc) = ^ iand < 9 (t)=<(*l. (65) 

Note that since a is the spread in quasimomentum, the 
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corresponding scaled variable is a 
troduced the scaled force [28] 

bF 



cj/Ul- We also in- 



(66) 



time (Bloch period) 



that compares the energy drop over a unit cell, bF, with 
the characteristic energy of the system, Eft. 

For the example shown in Figure 2 the ratios 
are t osc (0)/t b « 0.003 and T dccay (0)/T B « 0.004, 
while for the example shown in Figure 3 they are 
Tosc(0)/t b « 0.001 and T d ccay(0)/rs « 0.01. In both 
electron cases the ratios are small because F sa 0.002 
is also small, even though we employed a high electric 
field to provide the force. In contrast, for cold atoms 
in optical lattices it is easy to make these ratios bigger, 
while maintaining coherence over times of the order of a 
Bloch period. 

Consider, for example, the expectation value of 
the acceleration calculated for rubidium atoms pre- 
pared at the center of the lowest band N — of an 
optical lattice with b = 390 nm, s = 7 and accel- 
eration a,L = 24.2 m/ s 2 . The ratios (63) and (64) 
are now larger than in the electron case due to the 
increase in the force parameter to F sa 0.173. Here 
t OS c(0)/tb w 0.105 and r docay (0)/r B w 0.425, which 
correspond to r osc (0) sa 51.1/xs and Td OC ay(0) ~ 207/^s. 
Figure 4 shows how the atom's acceleration oscillates 
around the behaviour expected if the atom responded 
with the usual effective mass at all times. Similarly to 
the example shown in Figure 3, the acceleration goes to 
zero several times so we plot the acceleration instead of 
the time dependent effective mass (53). It is important 
to point out that expression (64) does not apply here 
strictly speaking because in this case the decay occurs 
as the center of the wavepacket has moved through a 
considerable portion of the Brillouin zone. Nevertheless, 
the result (64) still gives a rough idea of the initial decay 
time (see Figure 4). 

However, one of the most striking features of Figure 4 
is the revival of the oscillations as the wavepacket 
returns to k = completing a Bloch oscillation. This is 
not surprising when we consider expression (51), which 
predicts oscillations with frequencies and amplitudes 
periodic over the Brillouin zone. The revival is also 
present in the example shown in Figure 3 for the 
electron wavepacket when the acceleration is plotted 
for a full Bloch period. In the situation shown in 
Figure 2, however, we cannot use equation (51) to 
describe the behaviour of the electron wavepacket 
for an entire Bloch oscillation because the gap with 
respect to the band n = 3 (the closest band in the 
region near the edge of the Brillouin zone) becomes 
increasingly small as the wavepacket moves closer to the 
band edge (k = 1) and significant Zener tunneling occurs. 




40 80 120 160 
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FIG. 6. Expectation value of the velocity for a wavepacket 
(a = 0.2&l) initially at the center of the Brillouin zone in the 
band N = for a rubidium atom in an optical lattice with 
s = 7 and b — 390 nm, calculated from the first order approx- 
imation (red solid line), the full numerical solution (green 
dashed line) and the usual effective mass (blue dotted line). 
The acceleration of the lattice is 72.6 m/s 2 , which corresponds 
to F ~ 0.520. The inset shows the difference between the first 
order approximation and the usual effective mass prediction; 
this difference oscillates with a period approximately given by 
expression (63), r osc (0) « 0.315tb ~ 51.1/is. The units in the 
inset are the same as in the main plot. (Color online.) 



In the context of cold atoms, the wavevector ki, cor- 
responds to the wavevector of the optical lattice, and 
the energy Er, defined in (50), corresponds to the recoil 
energy, which is the energy gained (lost) by absorbing 
(emitting) one photon of the lattice. In (65) we used 
these quantities to define scaled variables for the wavevec- 
tor and the energy. We can also employ them to intro- 
duce additional scaled variables for position and time: 



x = k^x and t = — —t. 

n 



(67) 



Since the acceleration has units of force over mass, the 
corresponding scaled acceleration is [29] 



<a(*)> = ^-M*)) = S-("i''i>- «*> 



irEr- 



h 2 ir 3 



From (51) we have 



(~a(i))~F / dk\f N (k)\ 



rh* N (k) 



4 y2 t-^^ln PNnifypnNiK) COS7jV„(£, t) ) ■ (69) 
^ {£nN(K)y J 
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FIG. 7. Expectation value of the velocity for a wavepacket 
(a = 0.2fci,) initially at the center of the Brillouin zone in the 
band N = for a rubidium atom in an optical lattice with 
s = 13 and b = 390 nm, calculated from the first order ap- 
proximation (red solid line), the full numerical solution (green 
crosses) and the usual effective mass (blue dotted line). The 
acceleration of the lattice is 24.2 m/s 2 , which corresponds to 
F ~ 0.173. The inset shows the difference between the first 
order approximation and the usual effective mass prediction; 
this difference oscillates with a period approximately given by 
expression (63), r osc (0) ~ 0.0859tb ~ 41.8^is. The units in 
the inset are the same as in the main plot. (Color online.) 



Here we have introduced the additional scaled quantities: 



K — 


k- Ft, 


(70) 


In (a) = 


/f /jv(K) ' 


(71) 


Pnn' (&) = 


J-Pnn'(k), 


(72) 


7JVn(M) = 


-f Nn (K,t). 


(73) 



In many experiments that study the motion of cold 
atoms in optical lattices, the measured variable is the 
velocity of the wavepacket rather than its acceleration 
[3-8] . In Figure 5 we show a plot of the velocity for the 
same parameters as Figure 4. As before, we include a 
plot of the prediction based on the usual effective mass, 
which clearly shows a full Bloch oscillation. Note also 
the good agreement between the prediction based on 
the first order calculation discussed in section II and 
the full numerical solution, even for long times after 
the force was suddenly applied for both the acceleration 
(Figure 4) and the velocity (Figure 5). 

In the rest of this section we discuss the behaviour 
of the velocity in various situations realizable in experi- 
ments with cold atoms, where the natural scale for the 
velocity is given by the recoil velocity, 

fm 
mb 



Accordingly, we use the scaled velocity 

(«(*)) 



<«(*)) = 



VR 



(75) 



which results from integrating the dimensionless acceler- 
ation, 



(v(t))= / dt'(a(i')) 



VR 



(74) 



(76) 



Note that we assumed the wavepacket 's initial velocity 
is zero because, it starts at the center of the Brillouin 
zone in all the examples considered here. 

The deviations of the actual expectation value of the 
velocity from the usual effective mass approximation 
value are controlled by the amplitude of the oscillations 
of the scaled acceleration (69) and the recoil velocity 
(74). First, we fix the product of the mass and lattice 
constant to that used in the previous example (rubidium 
atoms with b = 390 nm) and explore the effect of 
changing (a(t)). The recoil velocity in this case is 
vr ps 5.89 mm/s. 

According to expression (69), (&(?)) for a given 
initial band N depends only on the width of the quasi- 
momentum distribution; the strength of the periodic 
potential, characterized by s; and the scaled force F. For 
fixed values of these three parameters the oscillations 
preserve their shape and are only rescaled when the 
physical parameters, such as the bare mass or the lattice 
constant, are changed. 

The width of the quasimomentum distribution and the 
reduced effective mass control how fast the oscillations 
decay according to (64). However, as pointed out 
before, the estimate (64) does not take into account the 
motion of the wavepacket through the Brillouin zone 
over one Bloch period since it was derived from the 
approximation (58), which misses completely the revival 
of the oscillations of the effective mass. The decay can 
be minimized or even removed if the oscillations of the 
effective mass do not have time to decay before the 
periodicity of the terms in (51) (over one Bloch period) 
returns the amplitude of the oscillations to the initial 
value. 

From (64) we would expect the decay time to become 
comparable to the Bloch period when increasing the 
scaled force F and the reduced effective mass m^(0). 
The first case is illustrated in Figure 6 where we show 
how the increase of F, keeping the other parameters 
fixed, eliminates the decay shown in Figure 5. The sec- 
ond case is shown in Figure 7 where, instead of changing 
F, we increase s in the parameters used for Figures 4 
and 5 so that the bands N = and n = 1 become 
flatter and the reduced effective mass increases. In both 
Figure 6 and Figure 7 the spread of the wavepacket is 
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FIG. 8. Expectation value of the velocity for a wavepacket 
(a — 0.004A;l, as has been achieved in [30]) initially at the 
center of the Brillouin zone in the band N — for a sodium 
atom in an optical lattice with s = 7 and b = 295 nm, cal- 
culated from the first order approximation (red solid line), 
the full numerical solution (green crosses) and the usual ef- 
fective mass (blue dotted line). The acceleration of the lat- 
tice is 800 m/s 2 , which corresponds to F ~ 0.173. The inset 
shows the difference between the first order approximation 
and the usual effective mass prediction; this difference oscil- 
lates with a period approximately given by expression (63), 
Tosc(0) ~ 0.105tb ~ 7.73^is. The units in the inset are the 
same as in the main plot. (Color online.) 



the same as for Figures 2-5. We can control the decay 
by changing the spread of the wavepacket as shown in 
Figure 8, where we plot our results for sodium atoms 
with F w 0.173 and s = 7 (as in Figures 4 and 5) for a 
very small spread of the initial wavepacket (a = 0.004). 
Notice that in this case, instead of a decay of the 
oscillations after half a Bloch period, there is a slight 
increase of the amplitude of the oscillations due to the 
higher sensitivity of the wavepacket to the changes of the 
momentum matrix elements and the energy differences 
at the center of the wavepacket as it moves through the 
Brillouin zone. Such small width in quasimomentum can 
be achieved experimentally in Bose Einstein condensates 
of sodium atoms [30] . 

The parameters F and s also control the amplitude 
of the oscillations of (o(f)) and (y(t)). From expression 
(69) it is clear that the amplitudes scale linearly with 
F (compare Figures 5 and 6). In the range of s values 
explored here (from s = 7 to s = 14), the effect of 
modifying s on the amplitude of the oscillations is much 
smaller than the effect of F (compare Figures 5, 6 
and 7). Thus, in an experimental setting, the oscilla- 
tions can be made more visible by increasing the force 
within the limits where Zener tunneling is not significant. 

The velocity deviations tend to be larger for lighter 
atoms and smaller lattice constants since the recoil 
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FIG. 9. Expectation value of the velocity of a wavepacket 
(a = O.Olfci) initially at the center of the Brillouin zone in 
the band N — 1 for a sodium atom in an optical lattice with 
s = 13 and b = 295 nm, calculated from the first order ap- 
proximation (red solid line), the full numerical solution (green 
dashed line) and the usual effective mass (blue dotted line). 
The acceleration of the lattice is 800 m/s 2 , which corresponds 
to F ~ 0.173. Parameters are as in [7]. The inset shows 
the difference between the first order approximation and the 
usual effective mass prediction. The oscillations shown in the 
inset start with a period approximately given by expression 
(63), r osc (0) w 0.176tb ~ 12.9/is; then they squeeze as the 
wavepacket moves through the edge k = kh of the Brillouin 
zone (where the gap between bands N = 1 and n — 2 in- 
creases); and finally they return to the starting period. The 
units in the inset are the same as in the main plot. (Color 
online.) 



velocity is larger. For instance, the expectation value of 
the velocity calculated for sodium atoms in an optical 
lattice with b = 295 nm for the same s, a and F as in 
Figure 5 is a simple rescaled version of the results in that 
figure according to the new recoil velocity of the system 
(vr ~ 29.4 mm/s). Some other examples of calculations 
for sodium atoms are shown in Figures 9 and 10 with pa- 
rameters previously used in experiments that investigate 
the acceleration of Bose Einstein condensates [6, 7]. The 
detection of the oscillations shown here imply resolving 
deviations from the usual Bloch oscillation that are small 
compared to the recoil velocity and probably within the 
uncertainty in the measurements in those references [31]. 

Figure 9 illustrates the case where the initial band is 
N = 1 instead of the ground state band used in all the 
other examples for cold atoms. An important difference 
with respect to the case with N = is the increase of the 
amplitude of the oscillations. Compare, for instance, the 
insets of Figures 7 and 9 in terms of the scaled velocity 
(75). Although Figure 7 corresponds to rubidium atoms, 
in scaled units it is equivalent to a plot of the velocity 
for sodium atoms with the same s = 13 and F ss 0.173 
but with the wavepacket starting in the band N = 
with spread a = 0.2. One of the reasons for this increase 
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is the smaller gap between bands N = 1 and n = 2 at 
/j = 0. Both figures predict almost no decay, but the 
reasons are different. For the situation in Figure 7 the 
cause is the small curvature of the band N = 0, which 
increases the reduced effective mass associated with 
N = and n = 1. In Figure 9, on the other hand, the 
cause is the small spread in quasimomentum and not the 
curvature of the bands; in fact, the larger magnitudes of 
the curvatures of N = 1 and n = 2 make the absolute 
value of the reduced effective mass associated with them 
smaller, and therefore the decay would be faster if the 
spread a were the same as for Figure 7. 

The last example, shown in Figure 10, corresponds to 
a situation where the wavepacket is again in the lowest 
band, N — 0, but under a force that is approximately 
two times larger than the one used in previous exam- 
ples with sodium atoms and produces an acceleration 
a = 1700 m/s 2 ; the strength of the potential is set to 
s = 14. Accordingly, the amplitude of the oscillations of 
the velocity is approximately doubled compared to the 
cases shown in Figure 8 and Figure 7 (for the latter the 
comparison is in terms of the scaled velocity). The first 
order approximation works very well due to the large 
gap between the bands N = and n = 1 for s = 14 and 
the deviations with respect to the full numerical shown 
in Figure 6 (s = 7) when increasing the force do not 
occur. In Figure 9 using s = 13 is not enough to prevent 
deviations with respect to the full numerical solution 
because the gap between the bands N = 1 and h = 2 
at k = is smaller than the energy difference between 
bands N = and n = 1. 

The validity of the first order approximation (69), 
relies on the populations in the bands n ^ N being small 
because, according to (46), there is no population of 
these bands to first order in A„jv(fc). In all the examples 
shown here the probabilities in the bands n ^ N do not 
exceed 2%, justifying the good agreement in most of the 
cases. The largest deviations occur for Figures 6 and 9, 
although in both cases the first order approximation 
describes correctly the overall behaviour. 

Figure 11 shows the populations of various bands as 
they change in time for the case presented in Figure 6. 
Notice that near the middle of the Bloch oscillation the 
population of the band n — 1 overshoots the population 
at the end of the Bloch cycle. This overshooting has 
been discussed before ([32], and references therein) when 
studying Zener tunneling for much stronger forces. The 
additional oscillations in the population observed in our 
case are related to the oscillations of the expectation 
value of the acceleration (and the velocity) around the 
prediction from the usual effective mass. Only a certain 
combination of amplitudes in the different bands yields 
a wavepacket behaving with the usual effective mass as 
discussed in Section II A (see expression (28)); thus, we 
would expect that the oscillations of the population are 
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FIG. 10. Expectation value of the velocity for a wavepacket 
(<r = O.Olfci) initially at the center of the Brillouin zone in 
the band N = for a sodium atom in an optical lattice 
with s = 14 and b = 295 nm , calculated from the first or- 
der approximation (red solid line), the full numerical solu- 
tion (green crosses) and the usual effective mass (blue dotted 
line). The acceleration of the lattice is 1700 m/s 2 , which cor- 
responds to F ~ 0.369. Parameters are as in [6]. The inset 
shows the difference between the first order approximation 
and the usual effective mass prediction; this difference oscil- 
lates with a period approximately given by expression (63), 
Tosc(0) ~ 0.176tb ~ 6.10/is. The units in the inset are the 
same as in the main plot. (Color online.) 



accompanied by oscillations in the dynamical response 
around the usual effective mass prediction. 

The connection between the oscillations in the popu- 
lation and the oscillations of the effective mass is more 
evident after comparing Figures 11 and 12. The latter 
shows the oscillations of the populations for the case 
discussed in Figure 9, which are faster as expected for 
faster oscillations of the velocity shown in Figure 9 com- 
pared to Figure 6. This relation between the oscillations 
of the population and the dynamical response is found 
in all the other examples, where the amplitudes are 
smaller than those shown in Figures 11 and 12. 

Figure 12 also explains the discrepancy between the 
full numerical solution and the first order approximation 
in Figure 9, near the end of the Bloch period. Com- 
pared to the example with rubidium atoms in Figure 6, 
where the main deviation occurs for the amplitude of 
the oscillations, the deviations in Figure 9 clearly show 
a superposition of faster oscillations. Such faster oscil- 
lations are a clear indication that the population of the 
neighbouring bands n — and n > 3 for the example 
in Figure 9 are more important than the population of 
the bands n > 2 for the example in Figure 6. Consistent 
with this observation, the probability in Figure 12 for the 
band n = 3 (< 0.30%) is larger than the probability in 
Figure 11 for the band n = 2 (< 0.07%). 
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FIG. 11. Probability of the rubidium atoms in the optical 
lattice described in Figure 6 to be in (a) the initial band N = 
(red solid line) and (b) the next two neighbouring bands: 
n = 1 (green dashed line) and n = 2 (blue dotted line). All the 
populations were calculated using the full numerical solution. 
Notice that the probability for the band n = 2 (< 0.07%) is 
significantly smaller than for the other two bands; we have 
not included the populations for higher bands since they are 
even smaller. The time is in units of the Bloch period. (Color 
online.) 




FIG. 12. Probability of the sodium atoms in the optical 
lattice described in Figure 9 to be in (a) the initial band 
N = 1 (dashed green line) and (b) the next three neighbour- 
ing bands: n = (red solid line), n = 2 (blue dotted line) 
and n — 3 (light blue slash-dotted line). All the populations 
were calculated using the full numerical solution. We have 
not included the populations for higher bands since they are 
smaller than the populations for the bands shown here. The 
time is in units of the Bloch period. (Color online.) 



V. CONCLUSION 

In summary, we have studied the effect of an exter- 
nal homogeneous force acting on a particle prepared 
initially in one band of a one-dimensional periodic 
lattice when the force is suddenly applied and remains 



constant afterwards. As predicted by Pfirsch and 
Spenke [2], the expectation value of the acceleration 
responds initially with the bare mass, and oscillates 
around the value predicted using the usual effective mass. 

Using the perturbation scheme of Adams and Wannier 
[10, 14] for a wavepacket located initially in one band 
only, we derived an expression for the expectation value 
of the acceleration that is valid over a full Bloch oscil- 
lation, provided that Zener tunneling is not significant. 
In this picture, the wavepacket responds with the usual 
effective mass as it acquires certain small components 
of Bloch functions from neighbouring bands, but due 
to the initial condition, which requires the wavepacket 
to respond with the bare mass, the expectation value 
of the acceleration (48) has various terms oscillating 
around the usual effective mass prediction (35) with 
different frequencies. For the cases considered here, 
the most important contribution in the sum over the 
different bands in (48) comes from the closest band 
to the initial one; accordingly, the frequency of the 
oscillations is governed by the energy difference between 
these two bands as the wavepacket moves through the 
Brillouin zone. The oscillations can decay because of 
the spread of the wavepacket in quasimomentum, but 
the periodicity of the momentum matrix elements and 
energy differences in (48) produces a revival of the 
oscillations as the wavepacket completes a full Bloch 
oscillation. 

We presented calculations for a toy model of a 
one-dimensional semiconductor illustrating the features 
described by Pfirsch and Spenke. The initial decay in 
this case is very fast (femtosecond timescale) and occurs 
while the wavepacket has moved over a small portion 
of the Brillouin zone, allowing us to write a simple 
expression for the envelope function that controls the 
decay of the oscillations according to the predictions by 
Pfirsch and Spenke (see expression 58). 

We also showed an analysis of the oscillations in a 
system of cold atoms in an optical lattice, where the 
timescales of the oscillations and the decay are much 
longer (of the order of microseconds) and comparable to 
a Bloch period. We analysed the effects of tuning the dif- 
ferent parameters, such as the force, the strength of the 
potential, the spread of the wavepacket in quasimomen- 
tum, the bare mass of the atom and the lattice constant. 
Since the velocity of the atoms is a good candidate for 
experimental measurements, we plotted its expectation 
value, showing how the deviations from the usual effec- 
tive mass prediction can be comparable to the oscillations 
calculated from the usual effective mass alone. Since the 
decohcrence in optical lattices is much smaller than in 
typical solid-state systems, it would be possible, in prin- 
ciple, to detect the oscillations during times of the order 
of a Bloch oscillation. In the case of optical lattices, the 
decay can be easily minimized or even suppressed when 
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the revival of the oscillations of the effective mass is faster 
than the decay due to the spread of the wavepacket in 
quasimomentum. This feature could be exploited to de- 
termine how much decoherence occurs during one Bloch 
oscillation. 
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Appendix: Wannier's procedure 

In this appendix we sketch the method developed by 
Wannier to decouple the bands to any order in the force 
F [14]. Equation (18) can be rewritten as 



(f) n (k,x) = W n (k)(f> n (k,x), 

(Al) 



for the modified Bloch states (14). Notice that equation 
(A.l) takes the form of an eigenvalue problem, but 
with the peculiarity that both the operator (acting on 
4> n (k,x)) and W n (k) depend on k. 

The parameter F introduced in (66) is appropriate 
to characterize how strong the external force is with 
respect to the lattice potential. For convenience, we 
use the other dimensionless variables defined in (65) 
and (67), and accordingly we introduce H = Ho/ Er 
and W n (k) = W n {k)/ Er. The Bloch states are kept 
unchanged so we can write ip n (k,x) = ip n (k,x) in the 
new variables. The same is assumed for the unitary 
transformation, U n 'n(k) = U n > n (k), and therefore 



n (k,x) = (p n (k,x). 

With these definitions (A.l) can be rewritten as 



d 

Ho - F ■ [ x + i — 
dk 



4> n {k,x) = W n (k)4> n (k,x) 



We attempt to solve this equation expressing 



and 



(A.2) 
(A.3) 

(A.4) 



as power series in F. It is assumed that the zeroth order 
corresponds to the usual Bloch states and band energies. 
Thus 



and 



tii%(k) = S n , n 



W^(k)=£ n (k). 



(A.5) 



(A.6) 



Replacing the expansions (A.3) and (A.4) in (A.2) and 
collecting terms with equal powers of F we can find a 
recurrent system of equations. To first order it is found 
that, in the original coordinates, the unitary transforma- 
tion U n , n {k) is given by (20) [10, 14]. Note that if F 
is small we expect the parameter |A n / n (fe)|, defined in 
(21) , to be small because the Lax connection is of the 
order of the lattice constant b, while the energy differ- 
ence is of the order of Er. The first order approximation 
for the energy W n (k) is simply the band energy renor- 
malized by the diagonal part of the Lax connection (see 
equations (11) and (22)). 
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